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Abstract. — With the help of exact ground states obtained by a polynomial algorithm 
we compute the domain wall energy A at zero-temperature for the bond-random and the 
site-random Ising spin glass model in two dimensions. We find that in both models the 
stability of the ferromagnetic and the spin glass order ceases to exist at a unique con- 
centration p c for the ferromagnetic bonds. In the vicinity of this critical point, the size 
and concentration dependency of the first and second moment of the domain wall energy 
are, for both models, described by a common finite size scaling form. Moreover, below 
this concentration the stiffness exponent turns out to be slightly negative 6s = —0.056(6) 
indicating the absence of any intermediate spin glass phase at non-zero temperature. 



Since Edwards and Anderson (EA) proposed the model for spin glasses, it has been 
discussed not only among researchers specialized for the subject, but in a rather wide com- 
munity of physicists working on random systems in general. Perhaps, this is largely because 
the EA model is the simplest possible model with short-ranged interaction for which we 
might expect a spin glass phase transition similar to that in the mean-field model, i.e. the 
Sherrington-Kirkpatrick model. However, a number of fundamental questions failed to be 
answered conclusively [ij. Even the very existence of the phase transition in three dimen- 
sions was questioned 0] recently. Whether or not the spin glass phase exists in the presence 
of a uniform magnetic field is even less clear || M. 

Another unanswered question is whether a spin glass phase (at non-zero temperature) can 
exist in two dimensions with an asymmetric bond distribution. It has been argued |5j, ||, 0] 
that an intermediate spin glass phase might be present in the p-T phase diagram between 
the ferromagnetic phase and the paramagnetic phase. In Fig. |], such a p-T phase diagram is 
shown including the proposed spin-glass transition line represented by the dash-dotted line. 
For the site-random model the evidence for the existence of a spin-glass phase seems to be 
even stronger than for the bond-random model || P| . 

On the other hand, in the case of the bond-random ±J model with p = 1/2, arguments 
for the absence of a spin glass phase in two dimensions were mainly based on results from 



Monte Carlo simulations O, [L2f and the estimates of the domain wall energy [13, 14 



The data from Monte Carlo simulations, however, are not available at very low temperature, 



e.g., below T = 0.4 J in Bhatt and Young's simulation |jlT[, which naturally made it difficult 
to exclude the possibility of the transition at a temperature smaller than 0.4J. Furthermore 
it is not clear, whether the "stiffness" exponent 6s is really negative because the data in 
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Figure 1: The schematic phase diagram with 

the previously proposed spin-glass transition 

f f) 

lines. Tc stands for the critical temperature 
of the Ising model on the square lattice. P, F, 
and SG stand for the paramagnetic, ferromag- 
netic and spin glass phases, respectively. 
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Cieplak and Banavar's paper jl4| clearly show a systematic positive curvature in a log- log 



plot of the domain wall energy versus system size for systems without vacancy. More recently 
results of Monte Carlo simulation at lower temperatures have been reported |]10| indicating 
a transition at T ~ 0.24J. In addition, 9s = was suggested |7] based on estimates of the 
domain-wall energy. These results are consistent with a finite temperature phase transition 
for which the low-temperature phase is only marginally or weakly ordered, meaning that the 
two-point spin-spin correlation function decreases algebraically as a function of the distance. 

The aim of the present letter is to reinvestigate this issue by studying the domain wall 
energy at zero temperature via the determination of exact ground states for large system sizes 
and huge sample numbers. This can be done very efficiently with the help of a polynomial 
algorithm described by Barahona et al. ||, which amounts to finding a minimum- weight- 
perfect-matching in a weighted graph with N = L 2 nodes and has computational complexity 
0(N 3 ). The model that we consider is the two-dimensional Ising spin glass with binary 
couplings denned by the Hamiltonian 

H = -^JijSiSj (1) 

where Si = ±1 are Ising spins, (ij) are nearest neighbor sites on an L x L-square lattice and 
the interactions strengths Jij are quenched random variables taking on one of the two values, 
+J and — J. We consider two different cases: In the bond-random model all interactions 
living on the bonds are independently distributed with a concentration p £ [0, 1] of ferro- 
magnetic bonds (J^ = +1). For the site-random model one first generates independently 
distributed random variables for all sites, q = ±1. The concentration of type-A-sites, i.e. 
those with e = +1, is c, the type-B-sites (e = —1) occur with probability 1 — c. Then, is 
set to be — J if and only if = e,- = — 1, and it is set to be +J, otherwise. In this case the 
ferromagnetic bond concentration is given by p = (2 — c)c. 

We calculate the domain wall energy A defined by A = E p — E a where E p and E a are 
the ground state energies with the periodic and the anti-periodic boundary conditions in 
the x-direction, respectively. Free boundary conditions are imposed in the ^-direction. Of 
crucial importance are the exponents p and 9g that characterize the system size dependence 
of the moments of the domain wall energy: 

[A] oc L p and [A 2 ] 1/2 oc L 9s (2) 

A positive value for p indicates the stability of a ferromagnetic ground state even in the 
presence of thermal fluctuations and thus the existence of the ferromagnetic long range order 



Figure 2: The scaling plot of [A] for the bond- 
random model, p^' = 0.896, <j)\ = 0.77 and 
ipi = —0.19 are assumed. The inset is the view 
focused on the region near p = pc . 
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at finite temperature |13| . On the other hand, a positive value for the stiffness exponent, 
9 S , with p being negative at the same time, still indicates the stability of the ground state, 
which possesses now long range order different from a ferromagnetic one. Thus a positive 
9s is interpreted as a sign for a spin glass phase at non-zero temperature. We define p^ 
and p^ as the critical concentrations of ferromagnetic bonds at which the asymptotic L 
dependences of [A] and [A 2 ] 1//2 , respectively, change from increasing to decreasing, i.e., the 
concentrations where a ferromagnetic phase and a spin glass phase, respectively, cease to 
exist at finite temperature. 

We computed [A] and [A 2 ] 1 / 2 for L = 4,6,8,12,16,24 and 32 at various values of p 
ranging from 0.50 up to 0.95. While the number of bond samples depends on L and p, it 
is 32768 for one of the most time consuming data points, such as the one for L = 32 and 
p = 0.5. We hypothesize the following finite size scaling form for [A] 

[A]L^ = h{{p-p^)L^). (3) 

The three parameters p^\ <pi and ipi have to be chosen such that a good data collapse 
for all data is obtained. To quantify the "goodness" of this fit, we used an appropriate 
cost function S(p c , 0, ip) || whose minimum value should be close to unity when the fit is 
statistically acceptable. When we use all the obtained data points, the best fit is achieved 
with S(p c ,<j),if}) ~ 2, which indicates that there is a non-negligible systematic error, i.e., a 
correction to scaling. Therefore, we have tried a similar analysis on a restricted set of data, 
omitting data which are presumably outside the asymptotic scaling regime, namely, data 
with p far from pW and data for small system sizes. For instance, the goodness of the fit 
can be significantly improved to S = 1.12 by omitting the data for p = 0.95, which yields 
the estimates 

= 0.896(1), 0i = 0.77(1), Vi = -0.19(2). (4) 

Considering the very small errors accompanying the data points, it is remarkable that 
all the data even including those for L = 4 can be expressed by the finite size scaling form 
(0). The resulting scaling plot is shown in Fig. |2|. We have confirmed that other choices 
of data points producing a value of S close to unity result in estimates of pW, 0i, and <pi 
consistent with the estimates quoted above. The value of pW is consistent with most of 
previous estimate such as 0.88(2) []T5|, 0.89(2) M and 0.89(1) |TB[ while inconsistent with 
0.885(1) i. 

In Fig. H [A 2 ] 1//2 is plotted against L. The lowest curve with crosses, which is almost 
straight with negative but very small slope, corresponds to p = 1/2. In other words, the 



Figure 3: The domain- wall energy [A 2 ] 1 / 2 of 
the bond-random model plotted against the sys- 
tem size L for various ferromagnetic-bond con- 
centration p. The inset is the view focused on 
the data points for p = 1/2. The straight line in 
the inset is obtained by the fitting to the data 
points excluding the two leftmost ones. 
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Figure 4: The scaling plot of [A 2 ] 1 / 2 of the 

(2) 

bond-random model. p c = 0.894, = 0.79 
and ip2 = —0.16 are assumed. 
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domain-wall energy decreases systematically but it does so very slowly. The method of 
least squares using all the data points yields 9s = —0.052(1) whereas all but the first two 
points (for L = 4 and 6) results in 9s = —0.060(2). Therefore, we quote here the value 
9 S = —0.056(6) as our estimate. 

Our results disagree with the suggestion 9s = —0.25 by Cieplak and Banavarfffl]. Con- 
sidering the size of actual reduction in [A 2 ] 1 / 2 as L grows from 4 to 32, we cannot rule out 
the possibility that the exact value of this exponent is 0, which means < lim^oofA] < oo. 
Such a scenario would be consistent with a suggestion by Ozeki0. In this case one has a 
marginal situation and we cannot decide whether the long-range order persists at a low but 
finite temperature based solely on a calculation of the stiffness exponent. We may, however, 
say that the low-temperature phase is only weakly ordered even if the phase transition takes 
place at a finite temperature. 

Similarly to the above-mentioned procedure employed for [A], we have tried a finite size 
scaling analysis for the data of [A 2 ] 1 / 2 , 

[A 2 ] 1 / 2 ^ 2 = f 2 {{p - pf^L^ 2 ). (5) 

Now we have to omit more data to get an acceptable fit with the value of S close to unity, 
indicating that the correction to scaling is larger for [A 2 ] 1 / 2 than for [A]. However, a good 
data collapse is obtained when we use only data for L > 12 and 0.85 < p < 0.91. The best 
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fit yields 

P M = 0.894(2), 2 = 0.79(6), ^2 = -0.16(4) (6) 

with S = 0.99. The resulting scaling plot is shown in Fig. |j. The present estimate of p^ 
is larger than but marginally consistent with all the previous estimates such as 0.86(2) fl~7f . 
0.85 and 0.870 @, while it is clearly inconsistent with 0.854(2) 0. 

It is remarkable that not only p^ but also 02 and ip2 agree with the corresponding values 
in (|4]) within the statistical errors. While the agreement in p c already suggests the absence of 
the intermediate phase, we consider the agreement in the critical indices as another evidence 
for the absence of the intermediate spin-glass phase, since it is hardly possible that the first 
and the second moment of A show the same critical behavior at different values of p c . 

We now focus on the site-random model. In Fig. |5|, [A 2 ] 1 / 2 is plotted as a function of 
the system size. A significant correction to scaling can be seen in Fig. |5|. The same remark 
applies to the first moment, [A]. We have performed a finite-size-scaling analysis similar to 
what has been done for the bond-random model. As for [A], in order to reduce the cost 
function, 5*, down to unity, the smallest system sizes L = 4 and L = 6 have to be excluded 
from the scaling plot. The data out of the range, 0.60 < c < 0.68, are also excluded in the 
quantitative estimation of c c and the indices 0i and ipi. As for [A 2 ] 1//2 , an even stronger 
correction to scaling is present as shown in Fig. |5], making an additional system size (L = 8) 
unavailable for the quantitative estimation of c c . The range of c from which the data are 
chosen is again 0.60 < c < 0.68. 

The critical concentration, c£\ and the critical indices, 0j and ipi, are defined in a similar 
fashion to @ and @, resulting in 

4 1 ) = 0.658(3), 0i = 0.78(2), ^ = -0.18(3) (7) 

4 2) = 0.661(4), 02 = 0.79(2), tfj 2 = -0.23(3) (8) 

Obviously cp = within the error bars, from which also we here conclude that no inter- 
mediate spin glass phase exists. In addition, again, the critical indices for [A] agree with 
those for [A 2 ] 1 / 2 within the statistical errors. These critical indices agree with those for the 
bond random model (([|) and @), implying that both models belong to the same universality 
class. 

To summarize, we have performed a systematic calculation of the domain wall energy at 
zero temperature with systems larger than previous calculations for both bond-random and 
site-random models. We observed a significant cross-over effect or a correction-to-scaling 
especially for the site random model, while no indication for a finite temperature spin glass 



phase could be detected. Some of the previous evidences for the positive stiffness exponents 
were based on systems smaller than those studied in the present paper and therefore may be 
attributed to this cross-over effect. The critical concentration for the ferromagnetic bonds 
and critical indices estimated from [A] agree with those from [A 2 ] 1//2 , again indicating the 
absence of an intermediate phase. Moreover, the critical indices for the site-random model 
agree with those for the bond-random model, which suggests that both models indeed show 
the same universal critical behavior and have qualitatively identical features away from p c . 
We have also seen that the domain wall energy seems almost independent of the system size 
below p c for the bond-random model at p — 1/2. Our result of a negative stiffness exponent 
9 S = —0.056(6) is statistically significant, although it is difficult to exclude the possibility of 
6s = because of the very small change in [A 2 ] 1 / 2 actually observed. 

Several other calculations are still in progress. We have calculated the spin glass suscep- 
tibility at very low temperature down to T = 0.05J, and again found no evidence for the 
phase transition at a finite temperature. We have also computed the domain wall energy 
with a replica boundary condition similar to the one used in Ozeki's work [rj|. These results 



will be published elsewhere f!9 |. As a future perspective it might be promising to perform 
a similar study of the two-dimensional Ising spin glass model with next-nearest neighbor 
interactions and for which indications of a spin glass transition have been found |pOf| . 

This work was mainly done while one of the authors (N.K.) was at the supercomputer 
center HLRZ c/o Forschungszentrum Jiilich, Germany. He would like to thank the HLRZ 
for its hospitality and financial support. 
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